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Table o 
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(2) 54 HISTORY ; Detailed Current Edit Histor 

(5) ° DECLARATIONS ; Declarative Part of ule 
(4) 1 MTHSHACOS = Standard H Floating Arc cosine 
(5) 18 MTHSHACOS_RB = Special HACOS routine 

(8) § $ MTHSHACOSD = Standard H Floating Arc cosine 
(7) 5 MTHSHACOSD_RB = Special HACOSD routine 
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: S FACILITY: MATH LIBRARY 
-+¢ 
: ¢ ABSTRACT: 
: ; MTHSHACOS is a function which returns the H ipo, point arc-cosine 
in radians of its H Floating point argument. The call is standard call- 
: ; by-reference, 
; ; MTHSHACOSD is a function which returns the H flection point arc-cosine 


in radians of its H Floating point argument. The call is standard call- 
; by-reference, 


; VERSION: 1 
: ; HISTORY: 
AUTHOR: 
John A, Wheeler, 20-Oct-1979: Version 1 
; MODIFIED BY: 


at ing Point Arc-cosine rou Sine SEP=1984 01: AX/VMS v04-00 P 
Deta led Current Edit History "$8 sats ri 9 32 38 ERTARTL. SR SRE ATH HHACOS.MAR; 1 — 3) 


-SBTTL HISTORY ; Detailed Current Edit History 


Edit History for Version 1 of MTH$HACOS 


001 = Adapted from MTH$GACOS version 1-002. JAW 20-Oct-1979 
- Added degree entry » pornts. RNH 22-MAR-1981 
Modified computation of 1 = x*2 to avoid loss of significance 
for arguments >=1/2 
- Eben e shared external references to G*. al Geaget=et 
minated s symbol lic short Literals. RNH 
bene MULH RO to MULH (SP)+, RO in MTHSHACOSD. RNH 20-Jan-82 
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Point Arc-cosine routine 
; Declarative Part of Modul 


~SBTTL DECLARATIONS 


: INCLUDE FILES: 
EXTERNAL SYMBOLS: 


OTSPARAMS .MAR 


GBL 

MTHSHATAN_R8 
MTHSHATAND R8 
MTHSHSQRT_R8 
MTHSSSIGNAL 
MTHSK_INVARGMAT 


; EQUATED SYMBOLS: 


none 


_MTHSCODE 


none 


1F 
“x898C 
“X92T1F é 
“x898(8469, 
*X68004007, “Xx 
*x0, a 


001, 
469, 
00 


> MACROS: 
> PSECT DECLARATIONS: 
.PSECT 
: OWN STORAGE: 
: CONSTANTS: 
H_PI_OVER_2: 
CONG 
* LONG 
H_PI: 
LONG 
“LONG 
H_90: 
LONG 
* LONG 
H_180: 
LONG 
“LONG 


“68004008, “x 
“x0, “x 


16-SEP- 
6-SEP- 


1384 91:37:18 aes Macro V04-00 Page (3) 
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; Declarative Part of Module 


Force undefineds to error 
Arctangent routine 
Arctangent routine 

Square root routine 

Math signal routine 

Error code 


PIC,SHR,LONG,EXE ,NOWRT 


Program section for math routines 


1.5707963267948966192313216916397514420 
3.1415926535897932384626433832795028841 
90 


180 
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i 
tandard H Floating Arc cos 6-SEP-1984 MTHRTL.SRCJMTHHACOS .MAR; 1 (4) 
-SBTTL MTHSHACOS = Standard H Floating Arc cosine 


oo 
oo 
“er 


p++ 
; FUNCTIONAL DESCRIPTION: 

; HACOS = H Floating point function 
; HACOS(X) is computed as: 


0, then HACOS(X 
1, then wht sty, 


: HATAN(SQRT(1=-X##2)/X 

then HA HATAN(SQRT((1-X) "(1+ 
6, then HA )/ 
178, then H «( 


rror. 


) = HATAN(SQRT(1-X##2 
X) = HATAN(SQRT((1=x) 


rPooo 
axe 
wu 


2 
< 
2 
/ 
< 


@tia 


CALLING SEQUENCE: 
hacos.wh.v = MTHSHACOS(x.rh.r) 
CALL MTHSHACOS(hacos.wh.r, x.rh.r) 
Because an H-floating result cannot be expressed in 64 bits, it is 
returned as the firs = argument, with the input parameter displaced 
to the second argument, in accordance with the Procedure Calling 
Standard. 

INPUT PARAMETERS: 


LONG = 4 3; Define Longuerg multiplier 
x = 2 * LONG ; Contents of x is the argument 


IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
hacos = 1 * LONG ; hacos is the result 


Bete Ge Ge Se Ge Ge Ge Ge Ge Ge Se Ge Ge Se Ge Se Ge Ge Ge Ge Ge Se Se 


00000004 
00000008 


00000004 
VALUE: H Floating arc-cosine of the argument 

IMPLICIT OUTPUTS: none 

COMPLETION CODES: none 

SIDE EFFECTS: 

Signals: MTH$S_INVARGMAT if iX! > 1 with reserved operand in RO/R1 (copied 

to sb mal mechanism vector CHFSL_MCH RO/R1 by ETessicna) . 


Associated message is: ‘INVALID ARGUMENT”. Result is reserved operand -0.0 
unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1. 


SOOOOCOCOOCOOSOCOCOOCOOOOOCOOCOCOCOOOOCOOCOCOOCO COO OOO OOO OOOO OOO OOOO OOO OOO 
a a a ee a a ee ec ee ee a ee ee ee ee ee dd 


NOTE: This procedure disables floating point underflow, enables integer 


eS SS . . 
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1-006 MTHSHACOS = Standard H Floating Arc cos 6=SEP=1984 11:24:2 MTHRTL.SRCJMTHHACOS .MAR; 1 (4) 
4 169 ; overflow. 
4 198 ; 
4 171 j;--- 
8 1g | 
41FC 04 174 ENTRY MTHSHACOS, “M<IV, R2, R3, R4, RS, Rb, R7, RB 
4 175 3; Standard call-by-reference entry 
§¢ 158 MTHSFLAG_JACKET ; Heath o chin’ toa’ dentet d i 
Be ¥ 3 Flag tha s is a jacket procedure in 
6D 00000000'GF 9E 4 MOVAB G*MTHSSJACKET_HND, (FP) 
4 ; set handler address to jacket 
nk 3 handler 
0049 178 ; Case of an error in routine 
0049 179 ; If an error, convert signal to user PC 
0049 180 ; And resignal 
Bote 181 MOVH @x(AP), RO ; RO/R3 = KX = @x (AP) : 
O4E 18@ BSBB  MTHSHACOS_R8 ; Call special HACOS routine 
0050 +=+18 MOVO RO, @hacos(AP) ; Store result in first argument 
Bpee 13 RET ; Return to caller 
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1-006 ATHSHACOS. RE - Special HACOS routine 6-SEP-1984 94:32: 38 EMTMRIL. § L.S oATHH HACOS.MAR; 1 ’ &, 
26 ] 7 -SBTTL MTHSHACOS_RB = Special HACOS routine | 
2° ! : ; Special HACOS = used by the standard routine and direct JSB call. 
056 191 ; CALLING SEQUENCE: 
056 136 ; save anything needed in RO:R8 ; 
056 193; MOVH 3 Input in RO/R3 
056 194 ; JSB AT ASHACO | 
056 195 ; return with mens fn RO/R3 
056 136 3 | 
bee 13 MTHSHACOS_R8 Special HACOS i 
3 i; Specia routine 
50 «BS 086 199 TST RO ; Fonpare X with 0 
06 12 5 9 BNEQ TEST_FOR_1.0 i; Branch if x = 
OSA 1 
005A 0¢ ; 
00SA 203; x= 0 
005A 04 ; 
Baan 05 
50 A2 AF 7DFD 005A 206 MOVO 4H_PI_OVER_2, RO ; RO/R3 = PI/2 
05 OOSF 07 RSB ; Return PI/2 if ix! = 0 
3088 08 
060 09 ; 
0060 10; 0< Ix! 
0060 11 3 
0060 \¢ 
0060 213 TEST_FOR_1.0: 
7E 50 70FD 0060 14 OVH RO, -(SP) 3; Save_X on stack 
50 8000 8F AA 0064 15 BICW #*x8000, RO 3 RO/R3 = {Xi 
08 50 71FD 0069 16 CMPH . ; Compare ‘Xi with 1.0 
18 006D 217 BGEQ  GEQ_T0_1.0 : Branch if !X! >= 1.0 
O06F 18 
O06F 19 ; 
O06F $0 :0< ixi < 1.0 
O06F 1; 
O06F ; 
50 OFAO 8F B61 Q06F CMPW #4000, RO 3; Check for possible loss of 
OF 14 0074 4 BGTR 1$ 3 significance 
7E 08 50 63FD 0076 5 SUBH3 RO, #1, -(SP) 3; (SP) = 1- xX 
50 08 60FD 0078 § ADDH #1, RO ; RO/R3 = 1+ X 
50 BE 64FD O07F MULH2  (SP)+, RO > RO/R3 = 1 = X42 
09 11 «0083 8 BRB 2$ ; Join main flow 
50 23 64FD 0085 9 1$ MULH RO, RO 3 RO/R3 = X#"2 
50 0 0 63FD 9 0 SUBH RO, #1, RO : orene = 1.0 - Xe#2 
00000000 ' GF 16 E 1 2$: JSB G*MTHSHSQRT_R8B : RO/RS = HSQRT(1-X*#22) 
6E 66FD 0094 ¢ DIVH2 (SP), RO ; RO/R3 = HSQRT(1-X*#2)/Xx 
00000000 ' GF 16 098 3 JSB G*MTHSHATAN_RB 3 RO/R3 = HATAN(HSQRT(1-X*#2)/X) 
05 11 +43 : BRB TEST_SIGN ; Branch to TEST_SIG 
OA § 3 
OA si<s ix 
OA 3 3 
OA 3 3 
OA 40 GEQ_T0_1.0: 
12 14 aa ie BGTR ERROR ; Branch to ERROR if {Xi > 1.0 
OA 43; 
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1-006 ATHSHACOS_RE - Special HACOS routine anes 94:32:38 MTHRTL.SRCJMTHHACOS.MAR; 1 ° (5) 
OA 64 ; Ix! = 1.0 
A 45 ; 
A 6 
50 7CFD A 4 CLRH RO ; RO/R3 = 0 
Qas 549 
oA 29 ; Test the sign of X in order to decide whether to add PI to the result 
ie 
8 A 535 TEST _SIGN: 
6E 73FD OOA 54 TSTH (SP) ; Test the sign of X 
06 18 Bpae 22 BGEQ RETURN ; Branch if xX > 
50 FF61 CF 60FD AA 6 apere H_PI, RO ; Add PI to RO/R3 if x < 0 
5E 10 0 4 57 RETURN: ADDL #T6, SP ; Clear stack 
05 008 58 RSB ; Return to caller 
0B4 59 
0B4 60 ; 
00B4 61; 1 < iki, error 
0084 $06 3 
00B4 6 
ae CO 0084 64 ERROR: ADDL2 #16, SP ; Clear stack 
DD 0087 65 PUSHL (SP ; Return PC from JSB routine 
OO'8F 9A 00B9 66 MOVZBL #MTHSK_INVARGMAT, -(SP) ; Condition value 
50 01 79 OOBD 267 ASHQ #15, #T, RO ; RO/R3 = result = reserved operand -0.0 
00C1 68 3 goes to signal mechanism vector 
43 69 ; CCHFSL_MCH_RO/R1) so error handler 
0C1 70 3 can modify the result. 
52 C 00C1 71 CLRQ R2 o aes 
Q00000000'GF 02 FB 00C3 ie CALLS #2, G*MTHSS$SIGNAL : Signal error and use real user's PC 
OOCA 27 : independent of CALL vs JS 
05 OOCA 274 RSB ; Return - RO restored from CHF$L_MCH_RO/R1 
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-SBTTL MTHSHACOSD = Standard H Floating Arc cosine 


+4 

; FUNCTIONAL DESCRIPTION: 

; HACOSD = H Floating point function 
; HACOSD(X) is computed as: 


If X = 0, then HACOSD(X) = 90. 

If X = 1. then HACOSD(X) = 0, 

If X = <1, then HACOSD(X) = 180, 

If 0 < xX £1, then HACOSD(X) = ATAN(SQRT(1-xX##2)/Xx). 

If 1/2. <x <'1 Shon HACOSD (Xx) : HATAND (SQRT((1-X) #(14X))/X). 

If -1/2 < x < 6, then HACOSD(X) = HATAND(SQRT(1-X##2)/X) + 180. 

I a x < 1/2, then HACOSD(X) = HATAND(SQRT((1-x)#(14X))/X) + 180. 
iMt, error. 


; CALLING SEQUENCE: 


hacosd.wh.v = MTHSHACOSD(x.rh.r) 

CALL MTHSHACOSD(hacosd.wh.r, x.rh.r) 

Because an H-floating result cannot be expressed in 64 bits, it is 
returned as the first argument, with the input parameter displaced 


to the second argument, in accordance with the Procedure Calling 
Standard. 


; INPUT PARAMETERS: 


LONG = 4 ; Define senquere multiplier 
x = 2 * LONG ; Contents of x is the argument 
IMPLICIT INPUTS: none 
OUTPUT PARAMETERS: 
HACOSD = 1 * LONG ; HACOSD is the result 
VALUE: H Floating arc-cosine of the argument 
IMPLICIT OUTPUTS: none 
COMPLETION CODES: none 


SIDE EFFECTS: 


Signals: MTH$_INVARGMAT if iX! > 1 with reserved operand ~ porns (co180ed 
to the signal mechanism vector ee MCH_RO/R1 by | 1B$S1G 

Associated message is: ‘INVALID ARGUMENT™. Resul 

unless a user supplied (or any) error handler changes CHFSL_MCH_RO/R1. 


NOTE: This procedure disables floating point underflow, enables integer 


is nee operand -0.0 
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cB 3; overflow 
cB 4; 
CB 5 ;--- 
cB § 

0cB 

4iFC O00CB 8 «ENTRY MTHSHACOSD, “M<IV, R2, R3, R4, RS, R6, R7, RB> 
CD 9 3 Standard call-by-reference entry 
OCD 40 ; Disable DV (and FU), enable Iv 
poco 41 MTHSFLAG_ JACKET ; Flag that this is a jacket procedure in 
6D 00000000'GF 9E Soc MOVAB G*MTHSSJACKET_HND, (FP) 

OD4 ; set handler address to jacket 

eee 3 handler 

Q0D4 8 ; Case of an error in routine 

00D4 4 ; If an error, convert signal to user PC 

00D4 44 : And dicesignal 

50 08 BC 70FD 00D4 45 MOVH ax (AP), RO ; RO/R = ax(A 
6 10 00D9 “6 BSBB MTHSHACOSD_RB ; Call special MACOSD routine 
04 BC 50 7DFD O00DB 4 mOVvO RO, @HACOSB (AP) ; Store result in first argument 
0 Bees $05 RET ; Return to caller 


50 

07 

50 FF36 CF 
7E_ 50 

50 8000 8F 
08 50 

31 

50 OFAO 8F 
OF 

7E 08 50 
50 008 

50 BE 

9 

50 0 

50.0 08) 50 
00000000 ' GF 
50. 6E 
00000000 * GF 
08 
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; Special HACOSD = used by the standard 


CALLING SEQUENCE: 
save anything needed in RO:R8 
MOVH RO 


— 10 
Point Arc-cosine routine 16-SEP-1984 91: 
- Special HACOSD routine -SEP=-1984 11: 


-SBTTL MTHSHACOSD_R8 = Special 


JSB MTHSHACOSD_RB 
return with result in RO/R3 


MTHSHACOSD_RB:: 
TSTw 


BNEQ D_TEST_FOR_1.0 


MOVO 4.90, RO 
RSB 


RO, -(SP) 

BICw  #*8000, RO 
RO, # 

BGEQ  D_GEQ_T0_1.0 


3s 0<¢ tnt < 1.0 


CMPW 94000, RO 


BGTR 

SUBH3 RO, #1, -(SP) 
ADDH RO 

MULH (SP)+, RO 

BRB $ 

MULH RO, RO 

SUBH3 RO. #1, RO 

JSB G*MTHSHSQRT_RB 
DIVH2 (SP), 

JSB G*MTHSHATAND_RB 
BRB D_TEST_SIGN 
x! 


1.0: 
Beat 10$ 
BRW ERROR 


32:48 th Macro V04-00 
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HACOSD routine 
routine and direct JSB call. 


; Input in RO/R3 


Special HACOSD routine 
Compare X with 0 
Branch if x = 0 


RO/P3 = 90 
Return 90 if ix! = 0 


; Save X on stack 

; RO/R3 = txt | 

3; Compare :Xi with 1.0 

; Branch if iX! >= 1.0 

; Check for possible loss of 
i; significance 

; (SP) = 1 =X 

; RO/R3 = 1+ X 

; RO/R3 = 1 = X*2 

; Join_main flow 

; RO/R3 = X22 

. RO/R = 1.0 ° Xe22 

; RO/R3 = HSQRT(1-X*e#2) 

3 RO/R3 = HSQRT(1-X*#2)/X 

3; RO/R3 = HATAND(HSQRT(1-X**2)/X) 
; Branch to D_TEST_SIGN 


Branch to ERROR if {X! > 1.0 


a 
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| 
131 4 : 
131 28 3 iki = 1.0 
131 410; 
131 411 
50 7CFD 0131 ei¢ 10$: CLRH RO ; RO/R3 = 0 
i d tia 
1 r $1? ; Test the sign of X in order to decide whether to add 180 to the result 
134 419 
134 ai3 D_TEST_SIGN: 
6 hed 13441 iB AL. (SP) i; Test ghe hie of x 
6 1 137 = 420 BGEQ yn TURN ; Branch 
50 FEF2 CF 60FD : ; : 1 5 RET Daas H_180, RO ; Add ‘$0 a *RO/R3 if x <0 
. a: eS 61 a : ce NipoL2 #16, SP 3; Clear stack 
05 014 424 RSB ; Return to caller 
014 425 
014 4 8 § 
014 4 3 1< Iki, error 
014 428 ; 
0143 429 
014 430 
014 431 
0143 3=432 .END 
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D_GEQ 10.1.0 18 S 1 
D-RETOR 15F R 1 
D_TEST_FOR_1.0 of R 1 
D-TEST_SIGN 154 R 1 
ERROR 4 R 1 
GEQ_10_1.0 AO R 1 
HACOS = 
waren? = 0 
4_180 0000 & 01 
as B28 & 01 
Nbloven. _ Soseuup 8 
Long~ = = s G00000b4 
MTIHSSJACKET_HND teeeeere 8 4 
MTH$$S1GNAL eereeere § 0 
MTHSHACOS 000 RG 01 
MTHSHACOSD QO0000CB RG 01 
MTHSHACOSD RB QO0000E1 RG 01 
MTHSHACOS RB 00000056 RG 01 
MTHSHATAND_R8& eeeeeeee =X 8 00 
MTHSHATAN_RB eeeeeses =X 00 
MTHSHSORT RB eeeeeeee x 8 00 
MTHSK_INVARGMAT eeeneeee x 00 
RETURN 00000080 R 01 
TEST_FOR 1.0 00000060 R 01 
TEST_SIGN OO0000A5 R 01 
y = 00000008 
ger es eee eee eee eee + 
! Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
eS, 0000000 ( :-} 00 ( 0.) NOPIC USR CON ABS LCL NOSHR NOEXE NORD NOWRT NOVEC BYTE 
_MTHSCODE 0000014 . Baber. Oe &- See PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
pa ee ee ee Oe ee + 
! Performance indicators ! 
tee ee eee eee ee ee eee mee ee ee 
Phase Page faults CPU Time Elapsed Time 
Initialization 31 0: 8:8 
Command processing 11 : 


Pass 1 8: 1 
Symbol table sort 0 :00: $ 
Pass 86 :00: 
Symbol table output 4 : _ 0 
Psect synopsis output 6 B: : Q 
33 ‘ 


Cross-reference output :00: 
Assembler run totals 0:00: 


-03 
01 
49 
-03 
01 


4 


wooo wuo—-eo 


Wo Wi COC 


0:00: 
3: 0: 
:00:1 
Tog, vort ing set Limit was 1050 pages. 
6438 bytes (15 pages) of virtual memory were used to buffer she intermediate code. 
There were 10 pages of symbol table space allocated to hold 27 non-local and 5 local symbols. 


492 source Lines were read in Pass 1, producing 14 object records in Pass 2. 


MTHSHACOS 
VAX=-11 Macro Run Statistics 


1 page of virtual memory was used to define 1 macro. 


; H Floating Point Arc-cosine ae A... ot 14 =138% 9: 3: g ee * Bacco raache man: 1 Page 13) 


bmw wenn ore tenn ene cere ewrereseoee > | 


! Macro Library Statistics : 


éqnqucen wee een eB eR Bee Dee BS + 


Macro Library name Macros defined 


_$255$DUA28: CSYSLIBISTARLET. MLB;2 0 
0 GETS were required to define 0 macros. 


There were no errors, warnings or information messages. 
MACRO/ENABLE=SUPPRESS1ON/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHACOS/OBJ=OBJ$:MTHHACOS MSRC$:MTHJACKET/UPDATE=(ENHS:MTHJACKET) #MSRC 
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